Jamming transition of a granular pile below the angle of repose 
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We study experimentally the relaxation towards mechanical equilibrium of a granular pile which 
has just experienced an avalanche and discuss it in the more general context of the granular jamming 
transition. Two coexisting dynamics are observed in the surface layer: a short time exponential 
decay consisting in rapid and independent moves of grains and intermittent bursts consisting in 
spatially correlated moves lasting for longer time. The competition of both dynamics results in 
long-lived intermittent transients, the total duration of which can late more than a thousand of 
seconds. We measure a two-time relaxation function, and relate it via a simple statistical model to 
a more usual two-time correlation function which exhibits strong similarities with auto-correlation 
functions found in aging systems. Localized perturbation experiments also allow us to test the pile 
surface layer receptivity. 



PACS numbers: 05.20.-y Classical statistical mechanics 
sandpiles; granular compaction 64.70.Pf Glass transitions 



45.70.-n Granular systems 45.70.Cc Static 



I. INTRODUCTION 

Dry, non-cohesive granular materials can flow like a 
liquid but can also sustain under the influence of grav- 
ity a finite angle of repose 9 r . Given the large num- 
ber of degrees of freedom involved in the dynamics, it is 
tempting to think in terms of phase transition: for a pile 
slope 8 > 9 r , surface layers flow downhill, whereas for 
< 9 r , these layers "freeze" and thepile behaves like a 
solid. As a matter of fact, Bak et al. proposed the idea 
that a granular pile "self-organizes" into a critical state at 
the angle of repose. However, a number of experimental 
works U3- Li LjJ- Lj- llJl have studied surface instability of a 
granular heap, namely avalanches, and have led to the 
conclusion that such avalanches do not behave in a crit- 
ical manner at all. 

In recent years, it has also been suggested that the 
jamming transition of granular materials could be the 
analog of a glass transition [?j,|8| . Further works on mean 
field glass models have developed this analogy and tried 
to unify concepts that had emerged in both fields, such 
as the dynamical temperature defined from fluctuation 
dissipation relations, and that related to Edwards' sta- 
tistical ensemble [1 El El III 03 

In this paper, we reconsider the avalanching pile prob- 
lem but, instead of focusing on avalanches statistics, we 
study the relaxation of the pile towards mechanical equi- 
librium following an avalanche. We first show that the 
dynamics exhibits non-trivial relaxations and is much 
slower than expected given the "microscopic" time-scale. 
We observe intermittent bursts which temporarily reac- 
tivate the pile activity and identify them with correlated 
movements in space and time. Perturbation experiments 
allow us to confirm the emergence of strong spatial corre- 
lations when increasing the pile slope. Finally we extract 
in the unperturbed case a two-time relaxation function, 
which can be shown to rescale in an aging like manner, 
including a 'reparametrization' of time . suggested by 



a simple model that we introduce in order to describe 
our experimental data. Moreover, this model allows to 
relate the above relaxation function to a two-time corre- 
lation function which shares strong similarities with auto- 
correlation functions appearing in aging systems. 

The experimental setup is described in section II 
whereas the observations are reported in section III. Sec- 
tion IV is then devoted to the presentation of the model. 
Finally, the experimental results are discussed in section 
V in the light of the model, emphasizing possible connec- 
tions with glassy dynamics. 



II. EXPERIMENTAL SETUP 

The experimental setup shown on figure ^ consists in 
a rotating drum of internal diameter D = 450mm and 
thickness S = 22mm, half-filled with steel beads of diam- 
eter d = 3 ± 0.025mm and density p = 7.8g.cm~ 3 . 




FIG. 1: Picture of the experimental setup. The inset displays 
a typical acquisition of the pile slope. 
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The friction coefficient and coefficient of restitution of 
the beads are respectively fi = 0.2 and e = 0.92. A DC 
motor drives the drum clockwise or anticlockwise at an- 
gular velocities as low as fi — 0.01°/s. One can fix the 
pile slope to any given value smaller than the repose an- 
gle by simply braking the drum. The pile is lightened 
by a continuous halogen lamp located far enough from 
the drum to prevent heating. Also the lab temperature 
is regulated at a constant value 20°C. The experimental 
area of interest is filmed either by a standard CCD cam- 
era (768 x 572 pixels at 25i/s) or by a fast one (480 x 210 
pixels at lOOOi/s) aligned along the axis of the drum. 
In both cases, the camera focus is set to have a spatial 
resolution of 1mm per pixel. Image acquisition and pro- 
cessing allows us to extract the pile slope, the surface 
roughness and the pixels where a displacement as small 
as 25/mi as occurred. 

The surface detection takes advantage of the gradient 
of the light intensity fluctuations which is maximum at 
the surface of the pile S(x). The averaged slope of the 
pile 9 is then deduced from a linear least square fitting of 
S(x) while the surface roughness r is given by the average 
deviation around the mean slope. Beads displacements 
are detected using an image difference method already 
described in In principle measuring the intensity 

difference between two consecutive images of the pile al- 
lows to identify the region where beads displacements 
have occurred. In practice, one needs to reduce the noise 
caused by lighting inhomogeneities, camera intrinsic vi- 
brations and acquisition noise. This is achieved by ac- 
quiring two consecutive sequences of n images: n — 25 
(respectively n = 200) for the standard (respectively the 
fast) camera. The two sequences are averaged resulting 
in two frames. Substracting them and binarizing them so 
as to remove any remaining noise leads to the detection 
of the area where displacements have occurred, n is cho- 
sen according to a compromise between the noise/signal 
ratio and the temporal resolution. The sensitivity of this 
detection method is calibrated by translating the camera 
with a micrometric translator: any displacement larger 
than 20/im is spotted unambiguously. The method takes 
advantage of the whole range of the acquisition dynamics 
to detect displacements much smaller than the pixel size. 
The counterpart is that only the existence of such small 
displacement is detected, but not its amplitude. 



III. EXPERIMENTAL RESULTS 

A. Repose and avalanche angles 

We first measure the repose and avalanche angles 9 r 
and 9 a of the granular pile. The drum is rotated at a 
very low velocity (fi = 0.01°/s). The pile slope then 
follows an intermittent regime of macroscopic avalanches 
HUE El- The pile slope 9{t) and the surface roughness 
r(t) are extracted every 5s from two experimental runs 
of 12 hours each. 




FIG. 2: (a): pile slope angle 9(t) (°) as a function of time, 
(b), resp.(c): histograms of 6 a to P , resp. Ostart- 



A typical sequence of the pile slope evolution 9(t) is dis- 
played on figure 121a) ■ In such a low rotation regime, 
9(t) increases linearly with time at the rate fi until an 
avalanche occurs at 9 s tart and the slope relaxes toward 
a stopping angle 9 stop - The pile roughness remains on 
the order of half a grain diameter and no correlations be- 
tween r(t) and 9{t) or their temporal derivatives emerge. 
Filtering the events for which the pile slope variation 
{9 s tart — 9 stop ) < 0.1° corresponds to a local slope rear- 
rangements implying only a few beads, we ended with 
a total number of a thousand of events identified as 
avalanches. The resulting histograms of 9 s t a rt and 9 stop 
- figure H£b),(c) - exhibit exponential tails, which could 
indicate non-trivial fluctuations in 9 start and 9 s t op - Still, 
in absence of a better definition, we choose to estimate 
the repose and avalanche angles 9 r and 9 a as the av- 
eraged values of respectively 9 s t op and 9 s t a rt leading to 
r = 19.2° ± 0.2° and 9 a = 19.7° ± 0.2°. The relatively 
small difference between 9 a and 9 r can be explained by 
the large restitution coefficient of our steel beads. 



B. Relaxation dynamics 

In order to investigate the relaxation process of the 
granular pile following an avalanche, the drum is rotated 
at constant angular velocity (l°/s) during a few minutes 
so that a number of successive avalanches have occurred. 
The drum is then stopped just after a given avalanche 
(9 = 9 stop ) and rapidly rotated backwards down to the 
chosen working angle 9 < 9 s t op , where it is braked firmly. 
The relaxation process is recorded with the standard 
CCD camera. Consecutive sequences of n = 25 images of 
the pile are taken every 15s. Beads displacements occur- 
ring in between two sequences are detected following the 
procedure described in section II. As a result one obtains 
a time series of binary images where black pixels indicate 
positions where displacements have occurred. 
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FIG. 3: The mobile volume fraction evolution 8A(t) for two 
different realizations at similar pile slope : (o) 9 = 15°; 
(+) 6 — 16.5°. Inset is the log-lin plot of the same data. 
Notice the exponential decay rate, which is identical in the 
monotonous case (+) and in the intermittent case (o). 
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FIG. 4: Displacements in the pile at time step (a) and (b) 
labelled on ngure[3] The dark pixels corresponds to positions 
where a displacement has occurred in the 15s preceeding the 
given time step. The grey overlay indicates the pixels among 
those, where displacements have occurred successively during 
30s following the given time step (see text for details). 



The "mobile volume fraction" of the pile SA(t) (i.e. the 
fraction of beads that have moved between the two acqui- 
sitions) is computed as the total number of black pixels 
normalized with the total number of pixels covering the 
pile. Figure |3| displays two typical records of 5A(t) as 
observed for two realizations of the pile relaxation with 
similar slopes 8 = 15° and 6 = 16.5°. During the very 
first time steps the relaxation process is identical in both 
records: the bulk of the pile relaxes rapidly from bot- 
tom to top on timescales of the order of the 15s delay 
between two consecutive images acquisitions. It involves 
isolated beads displacement on time-scale even shorter 
and of amplitude much smaller than the bead diameter. 
The relaxation process then slows down in a subsurface 
layer of thickness of the order of [10 — 20] beads diam- 
eters. In contrast with the bulk, the subsurface layer 



relaxation may differ significantly from one realization to 
another as it is the case for the two realizations presented 
here. In the following, we will focus on the relaxation of 
the subsurface layer. In one case, one observes a simple 
exponential decay of the subsurface layer activity with 
a characteristic timescale rj which in the present case is 
of the order of 200s. In the other case, the pile activ- 
ity evolution is more complex. One notices intermittent 
bursts which interrupt periods of exponential decay, with 
the same time-scale as in first case (see inset of fig [3J| . 
The competition between the exponential relaxation and 
the reactivation bursts results in an intermittent process 
occurring on much longer time-scales lasting up to 30 
minutes. 

A visual inspection reveals that the reactivation 
bursts correspond to collective motion of grain clusters 
whereas the exponential decay observed in between 
two bursts or when there is no burst at all involves 
individual beads displacements. This is illustrated on 
figure ^ which displays the inverse binary pictures of 
two time steps -(a) during a burst event, (b) during an 
exponential decay period. The mobile volume fraction 
illustrated by the relative number of dark pixels, 5 A is 
of the same order and their spatial distributions are very 
similar. On these pictures, the gray overlay indicates the 
areas, where the same pixels have successively record 
displacements during 45s (15s before and 30s after the 
given time step). In the case of the burst event, and in 
contrast with the exponential decay case, it shows very 
clearly, that those displacements which last in time are 
also correlated in space, in the form of a cluster of grains. 

We now investigate the influence of the pile slope on 
the relaxation dynamics by performing 185 realizations 
with o G [0° , 8 r ] . As shown on figure|Sl the probability of 
observing intermittent dynamics increases with the pile 
slope. No intermittent dynamics could be observed for 
o < 5° - most likely a statistical effect. 




FIG. 5: Probability of observing an intermittent relaxation 
dynamics as a function of the pile slope. The crosses show the 
realizations which have been performed and are set to for 
simple exponential decay and to 1 for intermittent dynamics. 
The vertical dash-dotted line indicates the angle of repose 
6 r = 19.2°. 
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Figure Efa) shows the averaged characteristic relax- 
ation time of the exponential decay t± as a function of the 
pile slope. It increases significantly from 50s for 9 — 0° to 
250s when 9 approaches the repose angle and can be fit- 
ted by a power-law in (9 r — 9)~ 1 . By contrast, the reacti- 
vation process time-scales do not exhibit significant vari- 
ations with 9. The distribution of time intervals between 
two successive bursts, shown on figure GJb), exhibits an 



exponential tail, P(T) 



-T/n 



where 77, = 100s is the 



typical waiting time between two successive bursts. In- 
terestingly one can notice that the pile slope for which 
the monotonous relaxation time rj, becomes of the order 
of the characteristic time interval between two bursts 17, 
coincides with an angle 9 m 10° for which the probability 
of observing bursts during the relaxation becomes signifi- 
cant (> 0.1). This indicates that bursts only occur when 
the relaxation process is still active in agreement with 
direct observations. 
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FIG. 6: (a) Averaged relaxation time for exponential decay 
realizations as a function of the pile slope and (b): Probability 
distribution of the time intervals between bursts in the case 
of intermittent dynamics. 



The relaxation being a non stationnary process, it is 
described more accurately by measuring a two-time re- 
laxation function. It is convenient experimentally to 
measure N(t w ,t) the averaged number of displacements 
per pixel between times t w and t w + t, given by: 



N(t w ,t) 



t w +t 



5A(u)du 



(1) 



where the brackets denote ensemble average over realiza- 
tions with characteristic time-scales of the same order. 
We compute N(t w ,t) for both types of relaxation. For 
the exponential decay case (figure Ufa)), we average over 
all realizations obtained for 9 < 5°. For the intermittent 
dynamics case (figure Efb)), we average over the inter- 
mittent realizations obtained for 9 > 13°. 
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FIG. 7: Nftujjt), averaged number of displacements per pixel 
between times t w and t w + t plotted vs. t increasing t w (top 
to bottom): (a) exponential decay case and 9 < 5°, (b) inter- 
mittent dynamics case and 9 > 13°. X: experimental data; 
continuous lines: 3 parameter fit for each set of curves (see 
text, section IV for details). 
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In the first case, N(t w ,t) rapidly saturates whatever the 
value oft w and remains smaller than one, which indicates 
that displacements occur in average less than once at any 
given position. This is in agreement with the picture of 
isolated and independent moves as described above for 
the exponential decay dynamics. In the second case, no 
saturation occurs. N(t w ,t) increases continuously and 
becomes much larger than one, which reflects the pos- 
sibility of several displacements at a given position. In 
neither cases, could we find an obvious rescaling of time. 
Figure \7\ also displays fitting curves of our experimental 
data. We will come back on these fits and their interpre- 
tation in section IV. 



C. Response of the pile to a localized perturbation 

In order to gather some indications on the state of 
the subsurface layer during the relaxation, we probe its 
response to a localized and instantaneous disturbance. 
During the pile relaxation, we let a bead fall on the pile 
surface from a height of 2cm, with zero initial velocity. 
The whole procedure is recorded by the lOOOi/s camera 
during Is covering the impact of the incident bead. The 
image difference method described in section II is ap- 
plied with two sequences of images acquired just before 
and after the impact in order to detect the displacements 
induced by the impact. We fixed the binarization thresh- 
old such that the very short delay of the measure ensures 
not to take into account displacements coming from the 
relaxation dynamics. Several realizations are performed 
for each pile slope and the reverse binary images obtained 
for each realization are averaged in order to obtain a gray- 
scale picture of the local probability of displacement. 

Two typical responses obtained for pile slopes respec- 
tively equal to 9 — 13.2° and 9 — 18.8° are shown on 
figure |HI a) and (b). Note that the depth of penetra- 
tion of the impacted area is essentially constant of the 
order of [10-20] beads diameters whatever the slope an- 
gle, the same thickness as the subsurface layer where the 
slowing down of the relaxation dynamics occurs. In con- 
trast, the lateral extension strongly increases with 9 and 
is larger uphill than downhill resulting in a asymmetri- 
cal response. The amplitude of the response is estimated 
by measuring the area of the impacted zone A{9) as a 
function of the pile slope. As shown on figure |S{c) the 
response increases strongly with 9 and can be fitted by a 
power-law in (9 r — 9)~ 5 . 

D. Summary 

We now summarize all experimental results. Just af- 
ter the occurrence of an avalanche the granular pile that 
we study relaxes very rapidly - less than 15s - in its 
bulk, but exhibits a much slower relaxation in a sub- 
surface layer of thickness [10 — 20] beads diameter. In 
this subsurface layer, exponential decay and reactivation 
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FIG. 8: Response in displacements of the pile to a bead im- 
pact. Grayscale pictures of the local probability of having a 
displacement for pile slopes to (a) 6 = 13.2° and (b) 6 = 18.8°; 
(the arrow indicates the localization of the bead impact, (c): 
Impacted area A{6) as a function of the slope angle. 



bursts compete resulting eventually in long-lived - more 
than 10 3 s - intermittent transients depending on the oc- 
currence of bursts during the relaxation process. At the 
bead scale, the exponential decay is composed of fast 
and small independent displacements, whereas the reac- 
tivation bursts are linked with correlated moves of beads 
clusters. The exponential decay time increases from 
50s to 250s like (6 r — 9) . The typical time interval sepa- 
rating two bursts Tfc ss 100s. Accordingly, the probability 
of having at least one burst becomes significant when 
becomes on the order of Tb, that is for 9 ss 10°. N(t w , t), 
the averaged number of displacements per pixel between 
times t w and t w +t behaves very differently depending on 
the dynamics which is considered - with or without reac- 
tivation events. Finally the instantaneous response of the 
pile to a given localized disturbance increases sharply in 
terms of spatial extension when the pile slope approaches 
9 r . These observations raise at least two questions of dif- 
ferent nature: 
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• How can one observe well defined time scales t± 
and Tb of the same order 10 2 s, given that the only 
characteri stic time that can be defined at the bead 
scale is y/d/g ~ 1 0~ 2 s? 

• Is there a simple model which allows a direct calcu- 
lation of N(t w , t) and proposes an interpretation of 
the relaxation dynamics of the pile? Could one find 
moreover a time rescaling for N(t w ,t), as it often 
happens in glassy dynamics? 

Answering the first question, which has to be consid- 
ered in the context of micro-mechanics, is out of reach 
of the present experimental setup. Still, recent numeri- 
cal studies 0, 0| reveal the existence of two types of 
constitutive structures, namely the strong forces chains 
network and critical contacts clusters. Whereas the char- 
acteristic size of the former remains constant, the size of 
the latter increases with 9. These results suggest an intu- 
itively appealing, although rather speculative parallel. A 
constant tj could be a fingerprint of the underlying role 
played by the strong force network. Conversely, the in- 
crease with 9 of the monotonous relaxation time t± might 
be associated with the growth of the critical clusters. 
Further numerical studies of these collective mechanisms 
are under progress. Also one should take into account 
the possible effect of aging at the contact scale. 

To deal with the second question, we introduce a sim- 
ple model in which particles - to be thought of as beads 
- can move under some specific conditions, such that the 
number of moves N(t w , t) can be computed analytically. 

IV. MODEL 



of their state - and sets them in an active state, with a 
probability v per unit time and per particle. If the chosen 
particle is already in a active state, this reactivation pro- 
cess has no effect. At this stage, the reactivation mecha- 
nism is considered not to involve any displacement of the 
bead, but rather a collective rearrangement of its envi- 
ronment. For this reason, reactivations are assumed not 
to contribute to the average number of moves N(t w ,t). 

Introducing the probability Pi(t) to be in state i at 
time t (i — . . . n), the fraction P m (t) of active particles 
is given by P m {t) — ^f =1 Pi(t). From this, we deduce 
the following evolution equation: 



— p = ~aP m + uP a = -aP m + v (1 - P m ) (2) 
at 

where the first term accounts for the total probability to 
relax towards the inactive state 0, and the second one 
corresponds to the reactivations starting from the inac- 
tive state. Let us now introduce the key ingredient of the 
model: thermal noise being irrelevant for granular media, 
it is reasonable to think that the reactivation process is 
driven by some mechanical noise generated by the relax- 
ation of beads in the vicinity of the particle considered 
(the range of 'interaction' might however be quite large, 
due to the presence of strong correlations, induced for 
instance by the force chains network) . This effect can be 
taken into account by imposing that the rate v of reinjec- 
tion depends, in a mean-field spirit, on the global number 
of active particles: v — v(P m ). We consider here only the 
simplest functional dependence: v = fJ>P m , where [i is a 
constant. Solving Eq. J2J using this assumption leads to: 



It has been argued above that the measured displace- 
ments were due to individual beads relaxation, which 
were from time to time reactivated by some collective 
rearrangements. In the following model, we further sim- 
plify this physical picture by considering that the du- 
ration of the collective rearrangement can be neglected, 
which amounts to say that the 'activity level' 5 A is in- 
stantaneously raised when a reactivation burst occurs. 
To implement this picture, we assume that particles can 
be in two distinct kinds of states, namely active states 
labeled from 1 to n and an inactive state labeled 0. In a 
active state, a particle can "move" spontaneously to any 
other state, with a rate (probability per unit time) a for 
transitions to the inactive state and a' for transitions to 
the other active states, resulting in a global rate of transi- 
tion starting from an active state 7 = a+(n—l)a'. These 
transitions are to be associated to physical displacements 
(or individual relaxations), and their average number is 
thus the analog of the experimentally measured quantity 
N(t w ,t). On the contrary, particles in the inactive state 
cannot evolve spontaneously, and thus do not contribute 
to N(t w ,t). In order to account for the collective re- 
activation process, we introduce an external mechanism 
which chooses particles at random - and independently 



Pm(t) 



I ding 



(1 dt 

where an auxiliary function g(t) has been introduced: 



(3) 



g(t) = (a - fi) + fiP m (0) 



1-e 



-(a-n)t 



(4) 



From these results, one can compute the average number 
N(t w ,t) of particle displacements (in the sense defined 
above) between times t w and t w + t: 



N(t w ,t) 



t m +t 



jP m (t)dt=l In (£%+^l 



The exact calculation of N(t w ,t) will allow us to fit our 
experimental data. In order to relate our results to the 
more general context of glassy dynamics, it is of inter- 
est to calculate in this model the correlation function 
C (t w , t) defined as the probability not to have changed 
state between t w and t w + t. This correlation decomposes 
into the sum of the probability not to move starting from 
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an active state at t w and the probability not to be reac- 
tivated given that the particle was inactive at t w . This 
leads to the following expression for C(t w ,t): 

C{t w ,t) = P m {t w )e-^ + [1 - P m (t w )\ ( 6 ) 

g{t w +t) 

We now check the validity of the model by fitting 
the experimentally measured N(t w ,t) with the expres- 
sion given in Eq. Whereas the above expressions 
depend on four parameters, a, 7, [A, the elementary tran- 
sition rates and P m (0), the initial fraction of active parti- 
cles, N(t w ,t) actually depends on three parameters only 
7p = 7P m (0), S = a — /i and \ip — /iP m (0). Extract- 
ing these parameters from the experimental data by a 
unique least-square fit of the full 2- variable measurement 
of N(t w ,t) for each type of relaxation (i.e. monotonous 
and intermittent) leads us to the fit displayed on fig- 
ures 0a), (b) and the experimental estimation of the pa- 
rameters written on the left side of table [H Figure 
explicitly demonstrates the quality of the rescaling ob- 
tained for N(t w ,t) as a function of g(t w + t)/g(t w ). 
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TABLE I: Parameters values extracted from the experimental 
data for each type of relaxation: (a) the exponential decay and 
(b) the intermittent dynamics. Left: direct estimation of the 
parameters by a unique least-square fit of N(t w ,t). Right: 
relaxation time scales 1/q and 1/5, and average number 77 of 
moves per beads (see text for details). 

We further estimate P m (0) by imposing the ratio 77 = — 
in the case of the simple exponential decay. This ratio is 
nothing but the averaged number of move per bead and 
can reasonably be set to 77 = 1 by considering that the 
exponential decay dynamics is precisely obtained when 
each bead simply relaxes once from the active state to 
the inactive one. The obtained value P m (0) = 0.54 is 
then imposed a priori to both dynamics. This allows us 
to determine the elementary timescales 1/a and 1/5 for 
both dynamics as well as r\ in the intermittent dynamics 
(see right of Table HJ) . These obtained values are in very 
good agreement with the direct observations: 1/a varies 
typically between 50s and 300s as given by the direct 
measurement of the exponential decay rate displayed on 
figEt a )- In the same way, 1/5 is perfectly representative 
of the total duration of the relaxation process, namely of 
the order of 1 /a in the case of the exponential decay and 
larger than 10 3 s when intermittent bursts occur. Finally, 
77 = 4.7 is a reasonable value for the average number of 
move per bead in the sense that it agrees well with the av- 
erage number of burst during an intermittent relaxation. 
Altogether, despite its simplicity, our model based on ele- 
mentary dynamical processes provides a good description 
of the relaxation dynamics. 
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FIG. 9: —N(t w ,t) as a function of g(t m + t)/g(t w ) for each 
type of relaxation: (a) the exponential decay and (b) the 
intermittent dynamics 

V. DISCUSSION 

Let us come back to the model properties and com- 
ment the analytical expression of the correlation function 
C(t w ,t), that we obtained. As mentioned above, this cor- 
relation is the sum of two contributions, a first one which 
decays exponentially on a characteristic time I/7, and a 
second one with a more complex form, which decays a 
priori on a longer time scale which depends on t w . Thus, 
one could expect C (t w , t) to exhibit the typical two-step 
relaxation form familiar to the glass community. How- 
ever, the first step of the relaxation is hard to evidence 
since its relative contribution is given by P(t w ) which 
vanishes as t w increases, so that we mainly focus on the 
second step of the relaxation. The corresponding term 
presents a particular scaling form, as it depends only on 
the rescaled function g(t w + t)/g{t w ) - this latter ratio 
also appearing in the expression of N(t w , t). Interest- 
ingly, this scaling form is precisely the gene ric aging form 
proposed by Cugliandolo and Kurchan [19j , which gener- 



alizes the simple aging scaling in (t w + t)/t w by introduc- 
ing a reparametrization of time g(t). A natural question 
in this context is to identify the characteristic time scale 
associated to the decay of h(t w ,t) = [g(t w + 1) / '.g^™)] -1 , 
for a given t w . This characteristic time to can be identi- 
fied with the inflexion point of h(t w , t) when plotted as a 
function of Int. Figure ITUI displays the relaxation of this 
aging-like term, for several t w and for two different sets 
of parameters. 

Several regimes have to be distinguished according to 
the value of t w . For t w <C 1/(5, one finds t = 1/fiPo + t w 
(which can be approximated by l/a + t w ), whereas for 
t w » 1/(5, t saturates to the value 1/8. On figure HUT a), 
the parameters are set to arbitrary values such that time 
scales are widely separated, so as to evidence the three 
regimes described above, namely to « l/a, to « t w and 
to w 1/(5. One clearly observes the simple aging behaviors 
(i.e. to ~ t w ) in the time window l/a <C t w <C 1/(5, the 
earlier and later time regimes being stationary. One may 
also notice the strong increase of the plateau value when 
t w reaches 1/8 due to the fact that g(t) saturates for large 
t instead of diverging. The scenario in the present model 
thus appears a bit different from the full aging case. 
Figure HOf b) displays the same function using the pa- 
rameters values given in table [I] in order to extrapolate 
to large times the experimental data. The results are 
rather similar, in particular concerning the behaviors of 
the plateau, even if the intermediate aging regime is much 
less obvious, because of the weaker time-scale separation 
between l/a and 1/6. 

This saturation has two important consequences: on 
one hand, aging must be interrupted beyond a certain 
time scale, and on the other hand the correlation function 
does not decay to zero, but rather to a plateau value de- 
pending on t w , which goes to 1 when t w — ► oo. Although 
this behaviors might seem surprising at first sight, its 
interpretation is clear: the global activity of the system 
vanishes at long times, so that decorrelation is harder and 
harder to achieve. Aging is interrupted in the present 
case because of a complete dynamical arrest, whereas 
in (thermal) glassy systems aging ceases if the system 
reaches Boltzmann equilibrium, thus allowing for further 
decorrelation. One can expect this behaviors to be typi- 
cal of athermal jamming systems which are not driven by 
an external noise, but rather by a self-generated one, and 
for which the dynamics dies at long times. In this frame- 
work, the vanishing of the dynamics and the resulting 
interruption of aging is then encoded in the saturation of 
the time reparametrization g(t). This may be considered 
as a general scenario for the long time dynamics of (non 
driven) athermal jamming systems, like spontaneously 
relaxing granular piles. 



VI. CONCLUSION 

In this paper, we have studied the behaviors of a gran- 
ular pile in a rotating drum after avalanches. We briefly 
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FIG. 10: h(t w ,t) as a function of \nt, for different t w for two 
set of parameter values: (a) l/a = 10s, 1/8 = 10 4 s and (b) 
l/a = 294s, 1/8 = 1250s. 

consider the statistics of the angle of the pile 9{t) in 
the avalanching regime, measuring in particular the his- 
tograms of the starting angle 9 sta rt and of the stopping 
angle 6 stop- We also characterize the effect of a localized 
perturbation, showing that the perturbed area strongly 
increases with the angle 9 at which the experiment is 
performed. 

Our main results concern the long time mechanical re- 
laxation of a granular pile just after a avalanche has oc- 
curred. Quantifying the 'activity' in the pile 8A(t) by 
comparing subsequent images, we show that two differ- 
ent kind of relaxations are found, namely monotonous 
exponential decays with a characteristic time of order 
10 2 s and intermittent relaxations for which spontaneous 
correlated moves, or activity bursts, generate long tran- 
sients lasting for more than 10 3 s. These long lived rear- 
rangements tend to be localized in a subsurface layer of 
thickness around 10 or 20 beads diameters. Interestingly, 
whereas the typical time r& elapsed between two bursts 
remains constant when varying the slope 9 of the pile, 
the exponential relaxation time T| depends strongly on 
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9, with a divergence of the form (6 r — d)^ 1 when 9 ap- 
proaches the angle of repose 9 r . Both kinds of relaxations 
can be found at least for 9 > 5°, and the probability to 
find an intermittent decay increases with 9. 

We have then introduced a two-time quantity N(t w ,t) 
measuring the average number of detected moves be- 
tween times t w and t w + t, which allows for a precise 
analysis of the intermittent relaxation process. Using 
the predictions of a simple stochastic model that we also 
introduce, with three fitting parameters, we find a very 
good rescaling of N(t w , t) for all values of t w , keeping the 



same set of parameter values. Theses values of the fitting 
parameters compare very well with what can be extracted 
directly from 5A(t). Moreover, computing a natural cor- 
relation function in the model in order to make contact 
with more familiar two-time quantity, we show that this 
correlation behaves essentially in an aging-like manner, 
as a function of g(t w +t)/g(t w ), where g(t) can be inter- 
preted as a reparametrization of time, which saturates 
for large times so as to account for the full dynamical 
arrest. 
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